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Abstract. 

We present a detailed analysis of the chemical abundances, ionization 
state, and origin of a partial Lyman limit system (iV(H I) ~ 10 16 cm~ 2 ) 
at low redshift (z = 0.167 towards PKS0405-1219). We also present the 
first estimate of the total baryons in extended gaseous envelopes around 
galaxies based on a simple photoionization model for the Lya absorbers. 
We find that extended gaseous envelopes can explain between 5-20% of 
the total baryons predicted by the big bang nucleosynthesis model. 
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Various studies have shown that the "forest" of Lya absorption lines observed 
in the spectra of background QSOs are a sensitive probe of the physical condi- 
tions of intervening gas (e.g. Rauch 1998). Comparisons of galaxies and QSO 
absorption systems along common lines of sight, therefore, allow us not only to 
determine the origin of QSO absorption systems but also to study properties 
of tenuous gas around galaxies at large distances. Over the past several years, 
studies of faint galaxies in QSO fields have demonstrated that galaxies of a wide 
range of luminosity and morphological type possess extended gaseous envelopes 
of radius ~ 180 kpc and that known galaxies of known gas cross section 
can explain most and perhaps all of Lya absorption systems of N(H I) > 10 
cm -2 (Lanzetta et al. 1995; Chen et al. 1998, 2001). A more detailed review 
of the relationship between Lya absorption systems and galaxies is presented 
by Lanzetta in these proceedings. Here we summarize the main reasons that 
support this conclusion as following: 

1. The galaxy-absorber cross-correlation function exhibits a strong ampli- 
tude on velocity scales of < 250 km s _1 and impact parameter scales of ^200 
kpc, and becomes diminishingly small at impact parameters beyond 200 h~ l kpc 
on all velocity scales (Lanzetta, Webb, & Barcons 1997). This indicates that 
Lya absorption systems strongly cluster around galaxies on dynamical scales of 
individual galaxies, rather than on scales of large filamentary structures. 
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2. There exists a strong anti-correlation between Lya absorption equivalent 
width and galaxy impact parameter. Namely, galaxies at larger impact parame- 
ters are found to be associated with weaker Lya absorption systems. Including 
the intrinsic luminosities of individual galaxies as an additional scaling factor 
further improves the anti-correlation between Lya absorption equivalent width 
and galaxy impact parameter, indicating that the extent of tenuous gas around 
galaxies scales with individual galaxy luminosity/mass (Chen et al. 1998, 2001). 

3. Almost all galaxies that are at Sv X]3000 km s _1 away from the background 
QSOs and at p^200 h~ x kpc produce associated Lya absorption lines, indicating 
that the covering factor of tenuous gas around galaxies is approximately unity 
(Lanzetta et al. 1995; Chen et al. 1998, 2001). On the contrary, almost all 
galaxies that are in the vicinities of the background QSOs ((5-u<3000 km s _1 and 
p < 200 hr 1 kpc) do not produce associated Lya absorption lines to within 
sensitive upper limits, consistent with the picture that tenuous gas around these 
galaxies are highly ionized due to an enhanced ionizing radiation intensity from 
the background QSOs (Pascarelle et al. 2001). 

By combining results of previous studies of extended Mg II, C IV, and Lya 
gas around galaxies (derived based on analyses of QSO absorption systems), we 
can establish a schematic picture of the structure of extended gas around galax- 
ies: A typical L* galaxy is surrounded by high column density gas of radius a 
few tens kpc (e.g. van Gorkom 1993), which, because of its high density, remains 
mostly neutral. The gas becomes more ionized as the gas density decreases with 
increasing galactic radius. Low-ionization species, such as singly ionized magne- 
sium Mgll, become the dominant observational signature out to rj 50 h~ l kpc 
(e.g. Bergeron & Boisse 1991). As the gas density continues to decline at larger 
galactic radii, the gas becomes still more highly ionized. High-ionization species, 
such as triply ionized carbon CIV, become the dominant observational signa- 
ture out to 100 h~ l kpc, at neutral hydrogen column densities N(R I) ps 10 16 
cm -2 (Chen, Lanzetta, & Webb 2001). The tenuous gas continues to extend to 
at least ps 180 h^ 1 kpc, at neutral hydrogen column densities at least as low 
as iV(H I) « 3 x 10 14 cm' 2 (Chen et al. 1998, 2001). This picture applies to 
galaxies of a wide range of luminosity and morphological type. 

There is, however, a lack of understanding of the physical properties of 
extended gas at large galactic distances, such as number density, ionization 
state, and metallicity, because of limited information regarding kinematics and 
metallicity of the gas. The primary difficulty arises in acquiring a high-quality 
QSO spectrum of wide UV spectral coverage. Here we present an analysis of the 
chemical abundances, ionization state, and origin of a partial Lyman limit system 
(LLS; iV(H I) « 10 16 cm" 2 ) at low redshift (z = 0.167 towards PKS0405-1219). 
Combining the results of our analysis for the z = 0.167 system with those for 
high-redshift Lya absorption systems obtained in the literature, we also present 
the first estimate of the total baryons in extended gaseous envelopes around 
galaxies based on a simple photoionization model. 

2. The Partial Lyman Limit System at z = 0.167 

The absorption system at z = 0.167 toward PKS0405-1219 {z ern = 0.5726) 
was first identified in an HST/FOS spectrum with a rest-frame Lya absorption 
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equivalent width of 0.65 A and a rest-frame CIV absorption equivalent width 
of 0.44 A. Subsequent echelle spectroscopy with HST/STIS further revealed ab- 
sorption features produced by H° (Ly/3), C + , N+, 0°, Si+, Si++, Si+ 3 , possibly 
Fe + and Fe ++ ; and most interestingly absorption doublets produced by N +4 
and +5 at the redshift of the absorption system. We estimated the physical 
properties of the absorption system based on both photoionization and colli- 
sional ionization models, using the relative abundance ratios between different 
transitions. Detailed analysis has been presented in Chen & Prochaska (2000). 
We summarize the results here. 

First, all the ions except N +4 and +5 showed consistent profile signatures 
that trace the partial LLS, while the NV and OVI doublets appeared to be 
broad with velocity centroids blue-shifted by ~ 30 km s _1 from the other ions. 
The differences in kinematic signatures strongly indicate that the N V and O VI 
doublets and the partial LLS traced by the other ions do not arise in the same 
regions (see Figure 1 in Chen & Prochaska 2000). 

Second, metals that are associated with the partial LLS showed a wide range 
of ionization state. We present the column density measurements in Table 1, 
which lists the ions, the corresponding rest-frame absorption wavelength Ao, and 
the estimated ionic column densities log TV" together with the associated errors 
in the first three columns. The comparable abundances between various Si ions 
indicate that the partial LLS is photoionized. 



Table 1. Measurements of Column Densities, 
Doppler Parameters, and Abundances 



Species 


Ao 


log N (cm- 2 ) a 


b (km s 1 ) 


[X/Hf 


HI... 


1025.72 


> 15.7 


34.7 ± 2.5 






1215.67 


> 15.7 






CII ... 


1036.79 


> 14.10 








1334.53 


14.27 ± 0.09 




-0.33 


Nil .. 


1083.99 


> 14.25 


11.6 ±1.0 


0.27 


NV .. 


1238.82 


13.84 ± 0.07 


50.2 ± 7.0 






1242.80 


13.91 ± 0.06 






OI ... 


1302.17 


13.68 ± 0.14 




0.25 


OVI .. 


1031.93 


14.67 ± 0.16 


72.6 ±4.9 






1037.62 


14.76 ± 0.07 






Sill .. 


1190.42 


13.22 ± 0.07 


9.9 ±0.9 






1193.29 


13.29 ± 0.05 








1260.42 


> 13.16 








1304.37 


13.40 ± 0.12 




-0.37 


Si III . 


1206.50 


> 13.33 




> -0.76 


Si IV . 


1393.76 


13.18 ± 0.05 








1402.77 


13.49 ± 0.07 




-0.45 



"Lower limits indicate that the lines are saturated. 
6 [X/H] is defined as log [N(X)/N(R)] - log [X/H] ( 



Third, comparisons of the relative abundances of the Si ions and predictions 
based on a series of photoionization models calculated using the CLOUDY soft- 
ware package (Ferland 1995) showed that the ionization parameter, which is the 
ratio of incident ionizing photons at 912 A to the total hydrogen number density, 
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OVI 



-6 -4 -2 

log U 

Figure 1 . Predicted relative abundances of various ions X for a plane- 
parallel slab of absorbing gas of solar metallicity based on a photoion- 
ization analysis using CLOUDY. The dashed lines indicate the plausi- 
ble range of the ionization parameter U determined using the observed 
abundance ratios between Si + and Si +3 for the system at z = 0.167; 
the dotted and dot-dashed lines indicated the lower and upper limits 
of U determined based on the ratios of Si ++ / Si + and Si + / Si ++ . 



U = 0912/ cnn, was constrained to be log U = —2.64 ± 0.07. In addition, the 
observed abundances of the N +4 and +5 ions cannot be explained by the same 
photoionization models (Figure 1). For an extreme upper limit log U = —2.2, 
both OVI and NV were predicted to be at least 0.4 dex less abundant than 
Si IV, while our measurements showed that OVI and NV are at least 0.4 dex 
more abundant than all the Si ions. This result, together with their kinematic 
signature, indicates that these ions arise in a different region and suggests that 
they may be collisionally ionized. 

Fourth, the temperature of the partial LLS was found to be T ~ 7.4 x 10 4 
K with a 1 a lower limit being T ~ 2.7 x 10 4 K and 1 a upper limit being 
T ~ 1.5 x 10 5 K, using the Doppler parameters b estimated from a Voigt profile 
analysis (column 4 of Table 1) for the unsaturated lines produced by Si + and 
N + . Because the cooling function is the most effective at T ~ 10 5 K (see e.g. 
Spitzer 1978), the temperature of the partial LLS is more likely to be a few 
times 10 4 K. 

Fifth, assuming that the more highly ionized component traced by the +5 
and N +4 ions is collisionally ionized, we estimated the temperature to be T = 
(2.6±0.1)xl0 5 Kby comparing the observed O VI and N V column density ratio 
with predictions based on a series of collisional ionization models presented in 
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Shapiro & Moore (1976). In order for the hot gas to remain at this temperature, 
we further estimated the mass of the galaxy or galaxy group that supports the 
extended hot gas to be ss 1.5 x 1O 12 M , assuming a simple isothermal profile of 
a half-mass radius r/j ~ 0.35 Mpc. By requiring a cooling time longer than the 
dynamical time, we estimated the number density of the absorbing gas should 
be < 3 x 1CT 5 cm" 3 . 

Finally, we estimated the elemental abundances of the partial LLS using the 
ionization fraction correction calculated from the CLOUDY model for an ion- 
ization parameter log U = —2.64. In addition, we also estimated the chemical 
abundance using the [0 1/ H I] column density ratio, which is a direct estimate of 
the oxygen elemental abundance in regions where the resonant charge-exchange 
reaction between oxygen and hydrogen becomes dominant. The results are 
shown in the last column of Table 1 for an H I column density log 7V(H I) = 16.5. 
It shows that the partial LLS would have 1/10 solar abundance determined for 
carbon and silicon, 1/5 solar abundance for nitrogen, and at least 1/2 solar 
for oxygen. Increasing (decreasing) iV(H I) by 0.5 dex decreases (increases) the 
abundance measurements by 0.5 dex. 

To summarize, the Lya absorption system at z = 0.167 is a chemically 
enriched partial LLS that exhibits a wide range of ionization state and is likely to 
be embedded in a more wide-spread, collisionally ionized hot gas that gives rise to 
the absorption features of +5 and N +4 . A summary of the physical parameters 
of the warm and hot gas that are associated with this system is presented in 
Table 2 of Chen & Prochaska (2000). Ground-based galaxy surveys (Spinrad et 
al. 1993; Ellingson & Yee 1994) of the field surrounding PKS0405-1219 have 
found an elliptical galaxy at z = 0.1667 and impact parameter p = 74.9 h~ l kpc 
with a rest-frame .ftT-band luminosity Lk = 1.20 Lk, and (2) a spiral galaxy at 
z = 0.1670 and at p = 62.8 h~ x kpc with Lk = 0.02 Lk,- By carefully examining 
deep optical and near-infrared images, we find that if there are unidentified 
dwarf galaxies at the absorber redshift and closer to the QSO line of sight, 
then they cannot be brighter than 0.02 Lk, • By applying the luminosity scaled 
anti-correlation between Lya absorption strength and galaxy impact parameter 
presented in Chen et al. (2001), we find that the elliptical galaxy alone can easily 
explain the amount of absorption. While it is not uncommon to find tenuous gas 
at large galactic distance (p ~ 100 h^ 1 kpc) around nearby elliptical galaxies 
(e.g. Fabbiano, Kim, & Trinchieri 1992) or in groups of galaxies (e.g. Mulchaey et 
al. 1996), it is surprising to find metal-enriched, high-column density gas at this 
large distance. The spectral characteristics of the elliptical galaxy exhibit signs 
of recent star formation (Spinrad et al. 1993), therefore the physical process that 
initiated the recent star formation might be responsible for transporting metals 
to large galactic distance. 



3. Total Baryons in Extended Gaseous Envelopes Around Galaxies 

Big bang nucleosynthesis (BBN) makes specific predictions for the light element 
abundances (D, 3 He, 4 He, and 7 Li) as a function of the nucleon to proton ratio rj 
(e.g. Schramm & Turner 1998). Abundance measurements of these elements de- 
termine r] and therefore reveal the mean baryon density £lb in the universe when 
combined with measurements of the cosmic background radiation. Measure- 
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ments of the deuterium abundance along three sightlines toward high-redshift 
QSOs have yielded a precise estimate of the baryon density of f2& h? = 0.02 (e.g. 
Buries, Nollett, & Turner 2001). Combined data from various cosmic microwave 
background experiments have also shown exciting agreement in the estimate of 

(e.g. Wang, Tegmark, & Zaldarriaga 2001). The consistent results in 
greatly increase the confidence in our current understanding of the formation of 
the universe. While all the baryons may be accountable by the Lya forest at 
high redshifts (Rauch et al. 1997), it has, however, been shown that the sum 
of all the baryons seen in different environment only makes up 1/3 of the total 
baryons predicted by BBN (Fukugita, Hogan, & Peebles 1998). Finding the 
missing baryons is therefore crucial for establishing a complete picture of the 
formation and evolution of the universe. 

Fukugita et al. (1998) attributed the "missing baryons" to the hot gas 
around galaxy groups with temperatures ~ 10 5 K, a temperature range that 
is out of reach by existing X-ray facilities but may be probed using the OVI 
absorption systems observed in the spectra of background QSOs (Mulchaey et 
al. 1996). Results of hydrodynamic simulations also suggested that 30%-40% 
of the total baryons are in ionized gas of temperatures between 10 5 and 10 K 
at low redshifts (Cen & Ostriker 1999; Dave et al. 2001). Recent analysis of 
intervening OVI absorption systems indeed showed evidence that supports a 
significant baryon reservoir in the hot gas probed by the O VI absorbers (Tripp, 
Savage, & Jenkins 2000). There is, however, a large uncertainty due to poor 
statistics in the line densities of the O VI absorption systems and unknown gas 
properties such as metallicity and ionization fraction of the gas. In addition, 
some fraction of the Lya forest have associated OVI absorption lines, a simple 
addition of the total baryon contributions from the two species will result in a 
double-counting error (see the contribution by Tripp in these proceedings). 

Because extended gaseous envelopes are a common and generic feature of 
galaxies of a wide range of luminosity and morphological type (Chen et al. 1998, 
2000), some fraction of the total baryons must reside in the gaseous envelopes of 
individual galaxies that have not been counted in previous measurements. Here 
we present the first estimate of the total baryons in tenuous gas around galaxies 
on the basis of known galaxies that are represented by a galaxy luminosity 
function and of known gaseous extent of these galaxies. Our analysis takes 
advantage of better known galaxy statistics and is based on hydrogen directly. 
It requires an accurate estimate of the ionization fraction of the tenuous gas, 
but it does not require knowledge of the metallicity of the gas. 

To determine the ionization fraction of an absorber requires high-resolution, 
high-sensitivity QSO spectroscopy to identify associated metal lines in the rest- 
frame UV spectral range. A handful of Lya absorption systems have so far been 
studied extensively for the ionization state of the absorbing gas (Bergeron et al. 
1994; Lopez et al. 1999; Prochaska 1999; Prochaska & Buries 1999; Rauch et 
al. 1999), in addition to the partial LLS at z = 0.167. All of these absorption 
systems are best explained by a simple photoionization model, in which the 
ionization fraction of the gas clouds only depends on the ionization parameter 
U. Adopting a best estimate of the ionizing flux intensity J912 for a gas cloud 
of known geometry, we can therefore infer the ionization fraction for a Lya 
absorber of a given neutral hydrogen column density. 
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Figure 2. Comparisons of estimated ionization fraction for seven Lya 
absorption systems and predictions based on a series of photoionization 
models. Points with circles indicate the two absorbers at z < 1. Curves 
indicate predicted ionization fraction as a function of neutral hydrogen 
column density iV(H I) for a range of ionizing flux intensity J912 at 912 
A and a range of neutral hydrogen clumping factor C. 



Figure 2 shows the predicted ionization fraction X = N(K + )/N(H) as a 
function of neutral hydrogen column density iV(H I) (curves) obtained using 
CLOUDY for a range of Jgi2- Estimates of a handful Lya absorption systems 
obtained from the literature are shown in closed points with the two z < 1 
systems highlighted in circles (the z = 0.167 system described in § 2 and the 
z = 0.79 system studied by Bergeron et al. 1994). We have also included a 

clumping factor C = J («Hi)/( n m) 2 ! where rini is neutral hydrogen number 
density, in the models to account for the dumpiness of the gas. Comparisons of 
various models show that a higher J912 would imply a lower C and vice versa 
for a system with known X and iV(H I). Given that the background ionizing 
radiation intensity evolves rapidly with redshift (e.g. Haardt & Madau 1996), we 
determine a best-fit range in J912 and C based only on the two z < 1 points. We 
find that it is unlikely to have J912 Z 10 -23 ergs sec" 1 cm" 2 Hz" 1 Sr" 1 and 10 
for these low-redshift absorbers. If we adopt the best estimated J912 obtained 
from the QSO proximity effect measurements (Pascarelle et al. 2001; Scott & 
Bechtold in these proceedings), 0.5-1 xlO -23 ergs sec -1 cm -2 Hz" 1 Sr" 1 , then 
we find a clumping factor of order unity. 

To estimate the total baryons in extended gaseous envelopes of galaxies, we 
first determine the total hydrogen number density profiles of individual galaxies. 
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log r (h _1 kpc) 

Figure 3. Hydrogen density profiles of a typical L* galaxy, assuming 
•^912 = 10~ 23 ergs sec -1 cm -2 Hz -1 Sr . Thick curves are for C = 1; 
thin curves are for C = 10. Solid curves indicate the estimated total 
hydrogen profile; dotted curves indicate the neutral hydrogen profiles. 
It shows that extended gas is mostly neutral at r < 15 h^ 1 kpc and 
therefore neutral hydrogen traces the total hydrogen gas profile. At 
r ^ 15 h^ 1 kpc, extended gas becomes highly ionized. While the neutral 
hydrogen gas profile appears to be as steep as nm cx r -6 - 5 , the total 
gas density profile is well represented by ran oc r~ 3 . 



Adopting the neutral hydrogen column density profile in Chen et al. (2001), 

/ n \ -5.56±0.42 / r \ 2.56±0.43 

7V H i(p,^) = 3.63xl0 21 ( -f— ) M- cm" 2 , (1) 



10 kpc/ XLb, 
we derive the underlying neutral hydrogen number density profile, according to 

nHl(rZfi) _ 1 f°° dN m (p,L B ) dp 

ttC J r dp ^/p 2 — r 2 

To derive a total hydrogen number density profile, we apply proper ionization 
corrections determined from the CLOUDY photoionization calculations for gas 
of a given nm- The results are shown in Figure 3 for extended gas around 
a typical L* galaxy. It shows that extended gaseous profiles may be well de- 
scribed by a power-law r a with a ~ —3, which is steeper than the one expected 
from an isothermal equilibrium model but agrees well with expectation from a 
dissipational formation scenario (Silk & Norman 1981). 
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3x 10 14 cm- 2 < JV(HI) < 10 19 cm- 2 
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Figure 4. Cumulative contribution of baryons in extended gaseous 
envelopes of galaxies of different luminosity. The calculations are for 
regions that correspond to 3 x 10 14 cm -2 ^ iV(H I) ^ 10 19 cnr 2 . 



Next, we calculate the cumulative contribution of baryons in extended 
gaseous envelopes from galaxies of different luminosity according to 

where /i is the mean atomic weight, toh is the proton mass, p c is the present-day 
critical mass density, and (j) is the galaxy luminosity function. The limits of the 
integral (n,r2) correspond to the column density interval of interest here with 
the lower limit iV(H I) = 3 x 10 14 cm~ 2 set by the sensitivity of the galaxy 
survey from which equation (1) was determined. 

Given the uncertainties in the galaxy luminosity function measurements, 
we calculate equation (3) using results from the CFRS (Lilly et al. 1995) and 
Autofib (Ellis et al. 1996) surveys. The results are shown in Figure 4. It shows 
that extended gaseous envelopes of luminous galaxies (Mb "Z Mb,) contribute 
most of the baryons, as indicated by the flattening of the curves toward fainter 
magnitudes beyond L*. The luminosity function of Lilly et al. has a higher 
normalization in galaxy number density and a shallower faint-end slope than 
that of Ellis et al., and hence gives a higher value of O gas - For a range of J912 
and C, we find that ~ 5-20% of the total baryons predicted by BBN are in 
extended gaseous envelopes of galaxies. 
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